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Abstract: We present a systematic study of the amplitude of the primordial per- 
turbation in curvaton models with self-interactions, treating both renormalizable and 
non-renormalizable interactions. In particular, we consider the possibility that the 
curvaton energy density is subdominant at the time of the curvaton decay. We find 
that large regions in the parameter space give rise to the observed amplitude of pri- 
mordial perturbation even for non-renormalizable curvaton potentials, for which the 
curvaton energy density dilutes fast. At the time of its decay, the curvaton energy 
density may typically be subdominant by a relative factor of 10~^ and still produce 
the observed perturbation. Field dynamics turns out to be highly non-trivial, and for 
non-renormalizable potentials and certain regions of the parameter space we observe 
a non-monotonous relation between the final curvature perturbation and the initial 
curvaton value. In those cases, the time evolution of the primordial perturbation 
also displays an oscillatory behaviour before the curvaton decay. 
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1. Introduction 

In the simplest curvaton mechanism |]l|, primordial perturbations originate from 
quantum fluctuations of a light scalar field which gives a negligible contribution to 
the total energy density during inflation. This field is called the curvaton a. Infla- 
tion is driven by another scalar, the inflaton 0, whose potential energy dominates the 
universe. After the end of inflation, the inflaton decays into radiation. If the infla- 
tionary scale is low enough, <^ 10~^y^, the density fluctuations of the radiation 
component are much below the observed amplitude Sp/p ~ 10~^ and the fluid is for 
practical purposes homogeneous. While the dominant radiation energy scales away 
as pr ~ a~^, the curvaton energy gets diluted at a slower rate, at least for quadratic 
curvaton potentials. The curvaton contribution to the total energy density therefore 
increases and the initially negligible curvaton perturbations get imprinted into metric 
fluctuations. The standard adiabatic hot big bang era is recovered when the curva- 
ton eventually decays and thermalizes with the existing radiation. The mechanism 
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can be seen as a conversion of initial isocurvature perturbations into adiabatic ones 
and, depending on the parameters of the model, is capable of generating all of the 
observed primordial perturbation. 

The scenario sketched above represents the simplest possible realization of the 
curvaton mechanism and a wide range of different variations of the idea have been 
studied in the literature. For example, the infiaton perturbations need not be neg- 
ligible 0, there could be several curvatons 0], the curvaton decay can result into 
residual isocurvature perturbations H and inflation could be driven by some other 
mechanism than slowly rolling scalars . Recently it has been pointed out that the 
curvaton can decay via a parametric resonance (see also [H) with potentially in- 
teresting observational consequences. It is also well known that the predictions of 
the curvaton model are quite sensitive to the form of the curvaton potential. In 
particular, even small deviations from the extensively studied quadratic potential 
can have a significant effect, at least when considering non-Gaussian effects ||^, p!0[] . 
However, thus far the predictions of a curvaton scenario with self-interactions have 
not been properly studied. This is the aim of the present paper. 

We explore in detail the curvaton scenario with the potential given by 

^=imV2 + A^-^, (1.1) 

where n is an even integer to keep the potential bounded from below and the in- 
teraction term is suppressed by a cut-off scale M. For non-renormalizable operators 
n > we set the cut-off scale to be the Planck scale, i.e. set M = Mp, and the 
coupling to unity, A = 1. For the renormalizable quartic case n = 0, we treat the 
coupling A as a free parameter. For the rest of the paper we set Planck mass to unity 
Mp = (87rG)-i/2 = 1. 

Our approach is purely phenomenological and we do not explicitly connect our 
discussion to any particle physics model in this work. Even without a concrete model 
in mind, the study of the potential ( |1.1|) is however well motivated by generic theo- 
retical arguments. Indeed, the curvaton should have interactions of some kind as it 
eventually must decay and produce Standard Model fields. The curvaton needs to be 
weakly interacting to keep the field light during inflation. This however only implies 
that the effective curvaton potential should be sufficiently fiat in the vicinity of the 
field expectation value during inflation but does not a priori require the interaction 
terms in ( |1 . ID to be negligible. Moreover, as the inflationary energy scale is relatively 
high, the field can be displaced far from the origin and therefore feel the presence 
of higher order terms in the potential. The interactions could arise either as pure 
curvaton self-interactions involving the field a alone or more generically as effective 
terms due to curvaton couplings to other (heavy) degrees of freedom that have been 



integrated out. An example of a possible physical setup which could lead to ( |1 . 1| ) is 



given by flat directions of supersymmetric models that have been suggested as cur- 



-2- 



vaton candidates [11]. This would lead to a potential of the form ( |1 . 1|) with typically 
a relatively large power for the non-renormalizable operator. 

We consider for simplicity a perturbative curvaton decay characterized by some 
decay width T and do not address the possibility of a non-perturbative decay [0]. 
Our analysis is numerical and covers both the regions of small and large interactions 
in ( |1.1| ) as compared to the quadratic part. For small interactions, the deviations 
from quadratic results appear mostly when considering beyond leading order prop- 
erties like non-gaussianties |P, |10] . The opposite case with large interactions is much 
less studied, although some work has been done |jl2|, |13|, and is therefore of pri- 
mary interest for the current work. When the interaction term dominates in (|1.1|) , 
the curvaton oscillations start in a non-quadratic potential and the curvaton energy 
density always decreases faster than for a quadratic case. For non-renormalizable 
interactions, the decrease is even faster than the red-shifting of the background radi- 
ation and the curvaton contribution to the total energy density is decreasing at the 
beginning of oscillations. Consequently, the amplification of the curvaton component 
is less efficient than for a quadratic model. For the same values for m and F, the 
curvaton typically ends up being more subdominant at the time of its decay than in 
the quadratic case. This we call the subdominant curvaton scenario. 

Despite the subdominance, the curvaton scenario can yield the correct amplitude 
of primordial perturbations as the relative curvaton perturbations (5(T*/cr* produced 
during inflation can be much larger than 10~^. For a quadratic model, it is well 
known that the curvaton should make up at least few per cents of the total energy 
density at the time of its decay, r = p^/pr ^ 0(10^^), in order not to generate 
too large non-Gaussianities This bound does not directly apply to the non- 

quadratic model ( |1 . 1|) since the dynamics is much more complicated. Although the 
subdominant curvaton scenario implies relatively large perturbations ba^ja^^ the 
higher order terms in the perturbative expansion of curvature perturbation can be 
accidentally suppressed P, [1^ . Although the study of non-Gaussianities is of great 
importance, stringent constraints on model building already come from the amplitude 
of CMB fluctuations. Therefore, in the present paper we focus on exploring the 
parameter space that produces the observed amphtude for the model ( p..l|) and leave 
the study of non-Gaussianities for a separate publication. As it will turn out, already 
this flrst step does place quite non-trivial constraints on the viable parameter space. 

We flnd that the observed 10~^ amplitude of primordial perturbations can be 
produced for all the interaction terms, n = 0,2,4,6, in (IJ.) which we consider. To 
illustrate the resulting constraints for the model parameters, we represent the corre- 
sponding value of r = po-/Pr at the decay time as a function of the inflationary scale 

and the curvaton energy during inflation r*. This describes how subdominant 
the curvaton ends up being at the time of its decay. For the special case of quartic 
self-interactions, n = 0, we also present a new analytical approach which can be 
used to estimate the curvature perturbation and compare this with our numerical 
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analysis. Our results imply that, as far as the amplitude of perturbations is con- 
cerned, the interacting model ( |1 . l]) is compatible with observations even for large 
interactions. This completes the first step in determining the observational status of 
the model ( |1.1|) and makes the study of the non-Gaussianities both well motivated 
and a necessary task to complete the work. 

The paper is organized as follows. In Sect. 2 we describe the subdominant 
curvaton model, the associated constraints, and discuss the qualitative features of 
the evolution of the curvaton energy density and the primordial perturbation. Sect. 
3 contains our numerical results. In Sect. 4 we present an analytical, approximate 
solution for the curvature perturbation in the case of a curvaton models which is a 
combination of a quartic and quadratic parts. Finally, in Sect. 5 we give a discussion 
of the results. 

2. The subdominant curvaton model 

2.1 Model constraints 

We consider a phenomenological curvaton potential ( |1 . 1| ) where m is the curvaton 
mass and the last term represents the leading interaction term. For n = we obtain 
marginally renormalizable four-point interaction. Larger exponents correspond to 
non-renormalizable operators and for these we set A = 1 and assume the associated 
cut-off scale is given by the Planck mass M = Mp = 1. We choose to consider even 
integer values for n to keep the potential bounded from below. The exponent values 
considered in detail in this paper are n = 0,2,4,6. For n > 6 there are no field 
oscillations in the non-quadratic part of the potential; as a consequence, are of no 
great interest. 

For the purposes of the present paper we assume for definiteness that the inflaton 
decays instantly to radiation so that p{4>end) = = pr. After inflation ends, the 
curvaton field must also eventually decay. We consider a perturbative curvaton decay 
characterized by a constant decay rate F, although the curvaton could also decay non- 
perturbatively by way of a parametric resonance 0. The potential ( |1.1|) then gives 
rise to a coupled set of equations of motion in a radiation dominated background: 

a+{3H + T)& + mV + A (n + 4) (t"+3 = 
pr = -4:Hpr + Fcr^ 
3H^ = p, + p„, (2.1) 

where pr is the radiation density and pa the curvaton density. We denote the relative 
curvaton energy density by 

r = ^ . (2.2) 

Pr 

The time evolution of r and therefore the rate of the amplification of the curvaton 
component depends on the order n of the interaction term in ( |1.1| ). In addition 
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to n, the curvaton mass m, and the decay width F, the solutions to the equations 
of motion depend on two initial conditions: the value of the inflationary scale if*, 
and the initial curvaton value a* which can be equivalently expressed in terms of 
the initial curvaton energy density r* ~ p^^/H"^. In principle, these are all free 
parameters albeit subject to the constraint coming from requirement of curvaton 
sub dominance during inflation r^, ^ 1 and masslessness V" <^ H^. Moreover, we 
assume that the curvature perturbations generated during inflation are negligible 
and that the evolution of if* is such that it gives rise to the correct spectral index 
for the curvaton perturbations. For example, a single field slow roll inflation with 
ii* ^ 10~^y/e^ and e* ~ 10~^ satisfies these conditions. 

An additional constraint comes from the time of the curvaton decay which should 
take place early enough so as not to spoil the success of the adiabatic hot big bang 
model. In particular, we require that the curvaton must produce only curvature 
and not (large) isocurvature perturbations. This corresponds to the situation where 
radiation and dark matter have both the same perturbation. Hence, we need to 
require that the curvaton decays before dark matter decouples from the primordial 
plasma. In other words, if dark matter decouples at the temperature Tdm, then 
F > iioM- For illustrative purposes, let us write F = g'^m with g 1, where g is 
some coupling constant, and assume Tdm ~ 0(10) GeV, as would be appropriate for 
LSP dark matter. Then we should require 

F > C(10)T2j^ ~ 10"^=^ ~ 10~^^ GeV, (2.3) 

or g > 10^^(TeV/m)^/^. In what follows, we will consider curvaton masses within 
the range m = 10^ . . . 10^^ GeV, so that for us the limit on the coupling constant g 
varies between g > 10^^^ . . . 10^^'', which from a particle physics point of view does 
not seem to be excessively restrictive. 

Although the initial value for the curvaton, and hence r*, is a free parameter 
constrained in our case by 

r* < 1 



one can obtain an estimate of its likely value based on a stochastic argument 
Consider the curvaton rolling down its potential during infiation. The comoving 
horizon shrinks and all the curvaton modes for which k > aH appear as stochastic 
noise in the dynamics, thus creating stochastic perturbations on scales k < aH as 
time passes. Note that, unlike for the infiaton in eternal infiation, for the curvaton 
there are no quantum "kicks" that would keep it indefinitely in the quantum regime; 
this is so because during infiation the background is fixed by the infiaton to be 
approximately deSitter and is almost independent of the subdominant curvaton. The 
one-point probability distribution P[(T{x),t] for the curvaton obeys a Fokker-Planck 
equation 

dP Id d'^P 



- 5 - 



The first term on the rhs arises from the classical dynamics while the second repre- 
sents the quantum noise. Here, by (t{x) we mean the curvaton smoothed on scales 
larger than the horizon and x denotes the position of one such superhorizon patch. 

At some point the comoving horizon will equal our observable universe and after 
that all the quantum noise of the curvaton will correspond to the perturbations we 
observe within our horizon. Therefore, the classical background value for the cur- 
vaton in our universe will be drawn by the solution to ( |2.4D at the time of horizon 
crossing for our observable universe. Of course this will depend on the initial condi- 
tion for the curvaton at the beginning of inflation, but, assuming that inflation lasts 
long enough, it is natural to consider the equilibrium solution of (|2.4|) . Thus, we 
have for the initial curvaton value a* 

Pequii = ^exp(^-||^l^(a.)^ (2.5) 

where H^, is the value of the expansion rate when our universe crosses the horizon. 
As a consequence, the characteristic initial value for the curvaton is given by the 
condition 

V{a) ~ ifVSTT^ (2.6) 

which is different from simply equating the classical force with the quantum fluctua- 
tions V ~ [|12[. Thus during inflation the subdominant curvaton has a probability 
distribution that is fairly flat. In Sect. 0, where we present our numerical results, this 
will not be considered as an additional constraint although we display the limiting 
value ( p^) in the Figures. 

Given these constraints, for each point in the parameter space and for each 
set of initial field values, we may then follow the evolution of the field and the 
expansion rate to compute the amplitude of the perturbation at the time when the 
curvaton finally decays. We use of the AiV- formalism 0, 0, where instead 
of following the evolution of background fields and perturbations around them, the 
evolution of superhorizon fluctuations is described by considering an ensemble of 
homogenous patches of universe with slightly different initial conditions. We compare 
the evolution of two such copies: one has the initial field value a*, while the other 
has 

cr+ = o"* + 5a^ , (2.7) 

where 5a^ = H^/2n. Since we are interested in the curvature perturbation on con- 
stant energy density hypersurfaces, these two patches should then be evolved to a 
fixed value of H, chosen so that for that value the curvaton has to a high accu- 
racy decayed completely in both patches. At this point the amount of e-foldings in 
both patches is compared, and the difference, AA^, gives the amplitude of curvature 
perturbation C, 

C = AA^ = A^((T+)-A^((T,) = Alna = (Ina)' 5a, + -(lna)" ^a,^ ^ . . . . (2.8) 



p 



2 



p 
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After the curvaton decay, the universe evolves adiabatically and the curvature per- 
turbation therefore remains constant C = on superhorizon scales. The amplitude of 
large scale primordial perturbations is fixed by CMB observations to be C — 4.8x 10~^ 



19, P 



2.2 Qualitative features of evolution 

Solving the full equations of motion analytically for arbitrary n is not possible, al- 
though for the special case n = some useful analytical approximations can be 



derived, as will be discussed in Sect. 4.1. Therefore, we will analyze the dynamics 



numerically. We tackle the full numerics in Sect. |], but let us first give a qualitative 
description of the dynamics using naive analytical estimates. The starting point here 
is the observation |^ that for a scalar field a oscillating in a monomial potential 



V ~ cr""'"^ with n < 6, the average energy density obeys the scaling law 

oc a"*'^ . (2.9) 



Let us therefore consider the limit where either the quadratic term or the interaction 
term cr""'"'^ dominates the potential ( |1 . 1|) . We find that there are three distinctive 
phases in the evolution of the curvaton. 

First, right after reheating, the curvaton is still effectively massless with V" <^ 
and the field value decreases very slowly. Once V" ~ if^, the curvaton becomes 
massive and starts to oscillate around the minimum of its potential. If the interaction 
term dominates at this stage, the oscillations begin in the a"'^^ part of the potential 
and the curvaton energy density scales as given by Eq. ( p.9| ). These oscillations 
continue, until |m^CT^ ~ Act"^'^, where a is the envelope of the oscillation. This 
marks the transition into curvaton oscillations in the quadratic part of the potential, 
whence the energy density scales as 

p^oca"^ (2.10) 

These three phases are clearly seen Fig. ||, which for a single choice of parameters 
displays the full numerical solution of evolution of r = pd Pr-, and then the piecewise 
scaling-law approximation, where we assume that when a single term dominates the 
potential, it scales according to Eq. ( p.9| ) until ~ F, and the curvaton decays 
inst ant aneously. 

The more subdominant the curvaton is during its decay, the less efficiently the 
initial perturbation in the curvaton field is transmitted to metric fiuctuations. How- 
ever, since the curvaton is very subdominant during infiation, i.e. its field value is 
small, the quantum fiuctuations produce relatively large perturbations when com- 
pared to the field value. Since the perturbation itself can be large, the curvaton does 
not need to be dominant at the decay. Instead, it is sufficient that its relative density 
grows enough from the given initial value. How much is "enough" will be discussed 
in detail in Sect. ^ 
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(a) Exact evolution. 
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(b) Piecewise scaling-law approximation. 



Figure 1: Comparison between the exact evolution and scaling-law approximation of the 
relative curvaton energy density. Here the parameters have been chosen to be n = 2, 
m = 10^^'^ and F = 10"^'^, and the inital conditions = 10^^^ and H^: = 10^". The three 
different phases of curvaton evolution can be clearly seen: First the energy density stays 
constant and hence r (x a^, then it starts to oscillate in the non-quadratic part and hence 
r DC a~^/^ until it enters the regime of quadratic oscillations where r oc a. Decay occurs 
when T = H. 



2.3 Breakdown of the scaling law and oscillation of the curvature pertur- 
bation 

From the qualitative discussion above, one would infer that the curvature perturba- 
tion varies smoothly and monotonously as the parameter values are changed. How- 



ever, this is not quite true as the naive use of the scaling law (|2.9|) does not capture 
all the relevant dynamical details. 

When using the scaling law (|2.9| ), one assumes the validity of two approximations: 
first, that a single term in the potential dominates completely; and second, that the 
oscillations are rapid when compared to the Hubble time H~^. Neither of these 
assumptions is necessarily valid. The presence of two different terms in the potential 



1.1) could be taken into account by generalizing the scaling law (p.9|) following the 



methods outlined in [^. This would not affect significantly the qualitative features 
of the results. The breakdown of the latter assumption of rapid oscillations can 
however have more significant effects. These can be seen in the numerical solutions 
of the equations of motion which exhibit complicated behaviour not qualitatively 
present in the simplified analytical approximation. To demonstrate this, in Fig. ^ 
we plot AN as a function of H for fixed values of if*, r* and F for two different 
curvaton mass values. For different masses the moment of transition from the non- 
renormalizable part of the potential to the quadratic part of the potential is different, 
and this affects strongly the final value of AA^ which is dictated by the duration of 
oscillations in the non-quadratic regime. 
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Figure 2: AA^ as a function of (time evolves from right to left). Both curves have 
-ff* = 10~^, n = 4, r = 10~^^ and r^: = lO"^'^, but masses are for the lower one m = 10~^^ 
and for the upper one 5-10"i^ Here we demonstrate the non-trivial evolution of AA^. Deep 
in the quadratic regime its oscillations become faster and faster so that AA^'s evolution 
is given by a scaling law. However, before the quadratic term in the potential starts to 
dominate, AA^ shows very complicated oscillatory behaviour. 

From Fig. ^it is clear that as the field value oscillates in time, so does AA^. In the 
non-quadratic regime AA^ oscillates with a large amplitude. If the transition to the 
quadratic regime is slow compared to the oscillations in the non-quadratic regime, the 
transition averages over several oscillations. As a consequence, the final value of AA^ 
will be a non-oscillatory function of the model parameters. However, if the oscillation 
frequency in the non-quadratic potential is slow, and the transition to the quadratic 
oscillations is rapid, then the phase of the non-quadratic oscillation affects the final 
value of AA^. If the parameters happen to be such that the transition to the quadratic 
regime occurs at a maximum of the oscillation, a relatively high value of AA^ freezes 
out. Similarly, if the transition occurs at a minimum of the oscillation cycle, the 
final value of AA^ will be much smaller. If the parameters governing the moment of 
transition, such as the curvaton mass m, are changed continuously, then the phase of 
the non-quadratic oscillation during the transition also changes continuously. In the 
space of the parameters this results in an oscillatory pattern in AA^. This behaviour 
can be understood by observing that the curvaton energy density at the beginning of 
oscillations can not be expressed in terms of an amplitude of the envelope alone but 
also depends on the phase of the oscillation, or equivalently on both the field a and 
its time derivative d, in a non-trivial way. In effect, these act as two independent 
dynamical degrees of freedom. If the transition from the interaction dominated part 
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to the quadratic region takes place at this stage, the initial variation of the curvaton 
value cr* can therefore translate in a non-trivial fashion into the final value of the 
curvature perturbation. 

In [1^ it was shown that for a potential V ~ cr"^'^ no oscillatory solutions exist if 
n > 6. This means that in the non-quadratic regime the curvaton merely decays and 
hence no oscillations in AA^ occur. This is consistent with the qualitative explanation 
of the oscillations discussed above. 



3. Numerical solutions 

3.1 Solving the equations of motion 

We have explored the parameter space numerically for subdominant curvaton models 
that produce the curvature perturbation ( ~ 10~^ using the AA^ formalism with the 
patches defined as in Eq. (p.7[). 

The solution of the full equations of motion ( |2.1| ) oscillates more and more rapidly 
as it enters the regime dominated by the quadratic term in the potential. As each 
oscillation cycle requires several steps, solving the full equations of motion becomes 
increasingly slow as the rate of oscillations increase. To resolve this problem, we use 
the fact that deep in the quadratic regime the oscillations are very fast and hence the 
scaling law approximation Eq. (|2.9|) becomes very good. Formally this corresponds 
to taking (d^) = p^. Then we may switch from numerically evolving the equations 
of motion ( p.l| ) to solving the energy density evolution equations 



p„ = -{3H + r) 

3H^ = p, + p^ . (3.1) 

Since the solution of these equations does not have oscillatory behaviour, the step 
size can be increased significantly without sacrificing accuracy. Thus we choose to 
evolve the full system until the quadratic term dominates the non-quadratic term 
by a very large factor, e.g. 10^ for n ^ 2. This value is chosen for convenience only 
and ideally should be as large as possible. This is the point where the numerical 
integration of the equations of motion becomes intolerably slow, and the smaller n 
in Eq. ( |1.1D is, the slower the numerics becomes. With n = 0, the system should 
be followed through a very large number of oscillation cycles. Thus for n = we 
have to stop intregration already when the quadratic-to-non-quadratic ratio is 10^ 
to prevent the time of integration to grow to several minutes. 

At this point, we switch to the evolution equations for the energy densities (|3.1| ) 
and follow the system until the curvaton has decayed, i.e. radiation dominates the 
energy density of the curvaton by a very large factor, which we fix to be 10^^. At 
this point we may compute the final value of AA^. 
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Figure 3: log^g ^decay as a function of (horizontal axis) and (vertical axis) for n = 0. 
The results are shown for a choice of coupling constant A = 10^"^. The white regions corre- 
spond to parameter values for which the amplitude of the curvature perturbation exceeds 
the observed value of ^ ~ 10~^, and are therefore immediately excluded. The character- 
istic value of the initial curvaton probability distribution V{a) ~ i/^/Svr^ is indicated by 
the dashed line with high probability below the line, while the area below the solid line 
corresponds to the condition for nearly scale invariant spectrum V" < 10~^H^. Below the 
dotted line non-quadratic terms in the potential are dynamically insignificant. 

We are interested in how subdominant the curvaton is when it decays. This 
subdominance is described by the ratio of curvaton-to-radiation densities at the onset 
of decay, defined as 

rdecay = ^(1 ~ H) . (3.2) 

Since the value of rdecay is not an exphcit parameter in the model, we calculate it 
by first solving the value for F which gives the final amplitude of the perturbations 
to be the observed magnitude, 10"^. rdecay is then calculated with these parameters. 
Hence the result of the numerical simulations is a value for rdecay and F for each 
parameter set {n, m, \, H^,r^}. 

3.2 Perturbation amplitude in parameter space 

In Figs. ||- ||we have plotted log^Q rdecay with different choices of n and m as a function 
of if* and r*. In addition to the value of rdecay we have plotted also the constraint 
for near scale invariance, V" = 10~^if^ (solid line). The dotted line corresponds 
to the equality of the quadratic and non-quadratic terms in the potential during 
inflation; below the line the non-quadratic term has no sizeable effect. We also show 
the position of the edge of the curvaton probability distribution as given by Eq. (|2.6| ) 
(dashed line); high probability is found below this line. 
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Figure 4: log^g ^decay cis a function of H^, (horizontal axis) and (vertical axis) for n = 2. 
The coupling constant is chosen to be A = 1. The white regions correspond to parameter 
values for which the amplitude of the curvature perturbation exceeds the observed value of 
C ~ 10"^, and are therefore immediately excluded. The characteristic value of the initial 
curvaton probability distribution V{a) ~ H^/Stt'^ is indicated by the dashed line with high 
probability below the line, while the area below the solid line corresponds to the condition 
for nearly scale invariant spectrum V" < W~^H^. Below the dotted line non-quadratic 
terms in the potential are dynamically insignificant. 

For n > we set A to unity, i.e. the cut-off scale to the Planck scale. We have, 
however, checked that decreasing A does not change the picture violently, but rather 
just moves the features seen in Figs. ^- §| slowly in the parameter space. For n = 0, 
A is treated as a free parameter up to the constraint A ^ 1 coming from the lightness 
of the curvaton in the interaction dominated regime. 

The colder colors reflect a more subdominant curvaton and the warmer more 
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(a) n = 4, TO = 10"* 



(b) n = 4, TO = 10 




10 10 10 

(c) n = 4, TO = 10"^^ 




10 10 

(d) n = 4, TO = 10-14 



Figure 5: Same as in fig ^ but with n = 4. 



dominant. The white color indicates that no matter how dominant the curvaton is, 
it cannot produce the observed amphtude of perturbations. The dominant curvaton 
scenario corresponds to the border between dark red and white, i.e. rdecay = 1- The 
numerical results demonstrate several features predicted by the simplified qualitative 



description discussed in sections 23 and 

The decay time is constrained from below by the requirement of no isocurvature 
perturbations, as discussed in Sect. ^TI|. Adopting the limit on F given in Eq. (|2.ci|), 
some of the parameter space is ruled out; this is the source of the cutoff at the lower 
left corner of the figures. 

Oscillations of the value of rdecay in the parameter space can be clearly seen. 
As discussed before, their source is the transition from the non-quadratic to the 
quadratic regime. Therefore, as expected, below the dotted line in the figures, where 
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(a) n = 6, m — 10 
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(b) n = 6, m = 10 



Figure 6: Same as in fig Q but with n = 6. 





10 10 

(a) logiQf decay 

Figure 7: log^Q r^ccay and log^Q F plotted as a function of H^: (horizontal axis) and r* 
(vertical axis) for m = 10~^^ and n = 4. Lines as in previous figures. 



the non-quadratic term never dominates, the oscillations of rdecay are absent. Like- 
wise, as there are no field oscillations in the n = 6 potential, there are no ensuing 
oscillations of rdecay in the parameter space, as can be seen in Fig. ^ Furthermore, 
below the dotted line the value of rdecay depends only on and not on H^. This is 
expected, since in the quadratic limit the relative curvaton perturbation is given by 



m 



5a^ _ H^/2tt _ 



Some of the parameter space of different curvaton models, as specified by n, and 
m, and by the inital conditions r^, and if*, is clearly ruled out. This is the white area 
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in the plots. Large parts of the parameter space are however allowed, and for these 
areas, we have calculated the degree of subdominance at the time of curvaton decay. 
The figures clearly demonstrate that a curvaton model is feasible even when it is 
subdominant by a factor of 10~^ during decay. However, one should note that a very 
subdominant curvaton usually implies a large inflationary scale, as is also evident 
in the figures. This might not be a problem but as the inflationary scale grows 
the curvature perturbations generated during inflation also grow and can become 
significant which needs to be taken into account in the analysis. For single field 
slow roll inflation, the inflaton perturbations are negligible if <^ 10~^y^. 

The results given by the numerical code described above were also checked with 
an independently developed program, to verify that the features seen are not numer- 
ical errors. 

4. Analytical results for n = 
4.1 Solving the equation of motion 

The numerical results discussed in the previous section can be contrasted with an 
analytical estimate in the special case of n = for which we the potential ( |1.1| ) reads 



For this particular potential, it is possible to find a reasonably accurate analytical 
approximation for the exact solution of curvaton equation of motion. To some extent, 
the analysis serves as a useful consistency check of our numerical results, but can 
also be of certain interest of its own as it allows in principle to treat arbitrary values 
of the coupling A analytically. 

In a radiation dominated background the equation of motion for a homogeneous 
curvaton with the potential ( [4.1|) can be written as 



where a = {t/t^Y^'^ is the scale factor, we have defined / = aa /a^, and the star denotes 
the end of inflation. The constants are defined as ci = A\al/Hl and C2 = rn^/Hl. 
To identify cr* with the curvaton value at t*, we should set the initial condition 
(j{t^) = (T*. However, as the field is light at early times and evolves slowly, we choose 
instead to formally set the initial conditions o"(t — ?■ 0) = a* and (T(t — )■ 0) = in the 
asymptotic past. This corresponds to 





(4.2) 




(4.3) 



and yields a(t^) 



a,(l - 0(ci)) 



cr*, assummg Ci > C2. 
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Exact solutions for Eq. ( [4.2|) can be found if one of the constants ci and C2 
vanishes. In the absence of the quadratic term, C2 = 0, the solution obeying the initial 
conditions ( [4.3[ ) is given by the elliptic sine fm=o{ci) = 2^/'^c^^^^sn{2~^/^a/^a, —1). 
This can be well approximated by the leading term of its trigonometric series 



-1/4^1/4 



Ci a 



(4.4) 



with b = it{2K{~1))^^ 1.1981. The approximative solution obviously satisfies 
a differential equation obtained by setting m = in Eq. (^.2| ) and making the 
replacement /(a)^ -> 6^c}''^2"^/^/(a). 

This observation motivates us to study a much simpler linear equation 



^+(c + x^)/(x)=0 



(4.5) 



instead of Eq. (|4.2|) . Here a yet another redefinition of variables has been made 

b^V2Xa, 



X = ac. 



1/4 
2 5 



m 



Instead of the initial conditions ( [4.3P we need to impose 

/(a;)|^=o = 0, f'{x)\:r=o = bc2^^^ 



(4.6) 



(4.7) 



to match with the desired exact solution of ( |4.2| ) in the limit m — )■ 0. By choosing 
the initial conditions ( [4.7| ) we implicitly restrict our analysis on the parameter values 
for which the quartic term is dominant or nearly dominant in ( [4.2| ) at early times, 
Ci ^ C2. Since Eq. ( [1.5[ ) yields the correct approximative solution in the quartic 
limit m — )■ and reduces to Eq. ( [4. 21) at late times x 1 when the quadratic 
term dominates the dynamics, we may expect it to yield a reasonable estimate for 
the asymptotic solutions of (|4.2| ). Comparison with numerical results confirms that 
indeed this is the illustrated in Fig. R|. As can be seen, the solutions of Eqs. 



c=l,c,=0 . 01 



c=10, Ct=0 . 01 





Figure 8: Comparison of the solutions for the curvaton field a = (T^f{a)/a obtained from 
the exact equation (E^) (green) and the approximative one (|4.5D (red). 



( [4.5| ) and ( [4.2p differ by a slowly varying phase but the amplitudes agree quite well. 
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This seems to be a generic outcome and not just a coincidence for the parameter 
values shown in Fig. For our current purposes the phase difference is irrelevant as 
it would affect the curvature perturbation at late times mt ^ 1 only through terms 
suppressed by m/H or r. 



4.2 Asymptotic solution 



The general solution for ( [4.5[) can be expressed in terms of parabolic cylinder func- 
tions which admit simple asymptotic representations in the limit of large argument 
[^ . Focusing on late times x = \/2mt ^ c, x ^ 1, Eq. (|4.5| ) yields the asymptotic 
result 



Cras(i) 



CTn 



(mt) 



sin ( mt 



TT 



6 



where 



(Jo 



|r (1 + ^1)1 



6 = |ln(2mt)-argfr^ + z^ 



(4.8) 

(4.9) 
(4.10) 



and the initial conditions ( [4.7|) have been used. Fig. § shows a comparison between 
Eq. (|4.8|) and the exact solutions of (|4.2|) for the same parameter choices as in Fig. 
P but for later time events x ^ c where the asymptotic solution is expected to 
be valid. Again it can be observed that the results are offset by a slowly varying 



c=l,c, =0 . 01 




c=10,c,=0 . 01 




Figure 9: Comparison of the exact solution (green) for the curvaton field a = a^f{a)/a 



obtained from Eq. (4^) and the asymptotic result (red) from Eq. ( [O 



phase difference, which manifests itself as the first term in Eq. ( [4.10| ), but the time 
evolution of the amplitude is well described by the asymptotic solution Eq. (^I8| ). 

The dependence of cXosc on the magnitude of the quartic term in Eq. ( [4.1|) is illus- 
trated in Fig |l^. For the interaction dominated case c ^ 1, Eq. ( [4.9|) yields Uosc — 



3/4 



a^c^^^'^2~^^% which according to ( [4.8|) leads to an estimate a ~ a^c~^^^{H/m] 
for the envelope of oscillations in the quadratic regime. This is in agreement with 
the result obtained using the simple scaling law (|2.9|) and choosing 4A(t^ 



as 
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^osc 




100 



Figure 10: The effective curvaton amplitude (|4.9D cJosc at the beginning of quadratic 
oscillations, c ~ 2\/Ao"*/m measures the relative strength of the interaction term at the 
end of inflation. 



the transition between quartic and quadratic regimes. In the limit c — )■ where the 
quadratic term dominates throughout the evolution, cTosc — > cr*6-\/7r2^^/^/r(3/4) and 
Eq. (|4^ ) reduces to an asymptotic representation of a Bessel function. Up to the 
constant b ~ 1.1981 related to the choice of initial conditions (4/7) engineered for 



the case c > 1, this coincides with the exact solution for the quadratic case obeying 
the initial conditions (14.31). 



4.3 Curvature perturbation 



Using the result ( |4.8| ), the curvaton energy density at late times mt 1, mt ^ c 
can be written as 



The contribution from the quartic part of the potential ( |4.1|) is negligible as equa- 
tions ( f4.8| ) and ( [4.9| ) yield Acr^/m^ < c^/(mt)^/^. Neglecting the non-gravitational 
interactions between the curvaton and the dominant radiation component, the total 
energy density reads 

p = pr*a~'^ + Pa , (4-12) 

where pr* — is the radiation energy density at the initial time t*. 

The curvature perturbation ( ^.8[ ) at some time t can be straightforwardly com- 
puted by varying ( [4.12| ) with respect to the initial curvaton value cr^, and keeping the 
final energy density p at t fixed. To leading order in r = Pa/Pr, this leads to the 
result 




oo 



C^'7;—Sa. + U^+{^] ]5al + ... = J2^Sa: , (4.13) 



where 
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72 2 /m\i/2 



and 

.-f-L(ij' . (4.15) 

The prime in ( |4.13| ) denotes differentiation with respect to a* and the amphtude of 
perturbations is estimated by 5cr* ~ H^j {2ti). In the sudden decay approximation, 
the curvaton is assumed to decay instantaneously at the time H = T where, up to 
factors of unity, V is given by the model dependent perturbative curvaton decay rate. 
The final value of the curvature perturbation under this assumption is obtained by 
setting H = Tm (|4l3|) . Here we assume that the universe evolves adiabatically after 
the curvaton decay and hence C = fl- 
it is also straightforward to derive a formally non-perturbative expression for 
the curvature perturbation (see also [Q]) which can be directly compared with our 
numerical results. By performing a finite shift a* — t- cr* + Act* in ( [4.12| ) but keeping 
the energy density p fixed, we arrive at a fourth order equation for the curvature 
perturbation ^ = In (a + Aa) — In a 



;i + r{a.)) S - 'y-*roscy-*-r^-*, - 1 = . (4.16) 



^osc(^*) 

4.4 Comparing with numerical estimates 

Equation ([4.16| ) can be readily solved using standard methods. Qualitative features 
of the solutions are illustrated in Fig. 11 below, which shows C('^*) for three different 
choices of the coupling A together with the corresponding numerical results. First, 
we can observe that the agreement between analytical and numerical results is very 
good, even surprisingly so. The mismatch seen for small field values occurs where 
the dynamics is dominated by the quadratic potential c <^ 1 and the choice of initial 
conditions for our analytical computation is no longer valid, as discussed above. In 
this regime, we expect an error proportional to 6 ~ 1.2 which seems to be well in 
line with the behaviour seen in Fig. |ll]. The bump around a^, ~ 10~® marks the 
regime where the quadratic and quartic terms are roughly equal initially, c ~ 1, and 
the transition from quartic to quadratic potential takes place soon after the onset of 
oscillations. For larger field values, oscillations in the quartic potential last for some 
while and the result can be understood using the simple scaling law ( |2.9| ). These 
features are generic although discussed here using the specific choice of parameters 
in Fig. |TT] as an example. 

Unlike for the non-renormalizable potentials, we do not detect any oscillatory 
behaviour of C in the regime c ~ 1 but instead only the bump. Indeed, using ( [4.16| ) 
we have checked that the first derivative of C, with respect to cr* is a positive definite 
function in the limit r <^ 1 and non-monotonous features will show up only in third 
and higher order derivatives. This differs qualitatively from the non-renormalizable 
case where C, oscillates as a function of a* and consequently even C' must change 
sign. Although we do not have an analytical estimate valid for non-renormalizable 
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10' 



10' 



Figure 11: The curvature perturbation ( for a choice of parameters n = 0, m = 10 , 
= W^'^ and T = 10"^®. The couphng values are A = 10"^ A = 2 x lO"'^ and A = 4x 10~'^ 
for red, green and blue curves respectively. The solid lines show numerical results and 
dashed lines analytical estimates. 

potentials, this agrees with general qualitative expectations as the transition into 
the quadratic regime takes place earlier due to the more rapid decrease of the field 
amplitude and is therefore more prone to leave an imprint on the perturbation. 

Finally, using ( [1.13| ) we can express the amplitude of Gaussian perturbations 
and the non-linearity parameter /nl = 5/6(lna)"/(lna)'^ to leading order in r and 
up to numerical factors in the form 

/m\i/2 1 /r\^^^ 

C ~ [y) ^'^""^ ' ~ ^ f -J X Uc) , (4.17) 

which coincide with the well known quadratic results apart from the functions 

fi = CToscC^osc/c^* and /2 = {aJaoscYil + o-qscC^osc/c^osc^)- These depend only on the 
dimensionless parameter c and their behaviour is depicted in Fig. [4.4| . It can be seen 
that the inclusion of a quartic term in the potential tends to decrease the perturbation 
amplitude and enhance the non-Gaussian effects, except for the slight decrease^ of 
/nl observed around c ~ 1. However, deviations from quadratic results are rather 
small as long as the quartic term is not vastly dominant over the quadratic part 

^Note that this differs from the result of Q based on a perturbative analysis. For n > 0, however, 
the sign change of /nl found in H is consistent with the oscillations of AN found in this paper. 
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Figure 12: /i and /2 plotted for different values of c ~ 2\/Acj*/m. fi and /2, defined in 



Eq. (|4.17D , describe the deviation of ( and /nl from the purely quadratic results. 



during inflation. For example, it can be seen from Fig. that the quadratic results 
get modified by at most an order of magnitude for c < 10, where c = 10 already 
corresponds to a clear dominance of the quartic term. In this regime the model 
is therefore compatible with observations under much of the same conditions as a 
quadratic model with the same parameters (a*, if*, m, F), although one also needs to 
take into account the different consistency constraints coming from the masslessness 
and subdominance of curvaton during inflation. If the interaction becomes more 
dominant, c > 10 or so. Fig. [4.4| implies that the main constraints come from the 
observational upper bound |/nl| ^ 100 for the bispectrum [^. This translates into a 



lower bound for the curvaton subdominance r, which will be more stringent than the 
quadratic result r > 0.01 [Q. These limits from non-gaussianity need to be checked 



more carefully, and calculated also for the cases > p3 



5. Conclusions 

We have studied the evolution of the amplitude of primordial perturbation in a class 
of curvaton models where the potential contains, in addition to the usual quadratic 
term, either a quartic or some non-renormalizable term. For the calculation of the 
perturbation amphtude, these turn out to be dynamically important for a large 
part of the parameter space. We also focused on the possibility that the curvaton 
energy density is subdominant at the time of its decay. To some extent, this a 
natural consequence of non-renormalizable curvaton potential since then the relative 
curvaton energy density r gets damped after inflation when the field is oscillating in 
the non-quadratic part of its potential. Depending on the decay time, this stage is 
followed by oscillations in a quadratic potential whence r starts to increase again, 
but typically remains subdominant for realistic decay widths. 

Although the amplification of curvaton energy is less efficient than for a quadratic 
model, the interacting model may still generate the observed amplitude of primordial 
perturbation since the initial curvaton perturbation (5o"*/cr* can be much larger than 
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10~^. This translates into constraints for the model parameters: the decay width F, 
the inflationary scale if*, the bare curvaton mass m and its initial energy density 
r*. In the present paper we have performed a systematic study of the available 
parameter space for different types of non-quadratic potentials. We have shown that 
the correct amplitude of perturbations can indeed be achieved for large regions of 
the parameter space even if the interaction term dominates over the quadratic part 
and the curvaton would be subdominant by a factor of r ~ 10~'^ or so at the time of 
its decay. 

Moreover, we found that the final relative curvaton energy density rdecay required 
to give rise to the observed amplitude depends on the model parameters in a highly 
complicated way and, for non-renormalizable potentials, shows oscillatory behaviour 
when the non-quadratic part dominates initially. The same oscillatory features show 
up when considering the curvature perturbation ^ as a function of the curvaton value 
(T* during inflation. Again, this happens for a large portion of the parameter space. 
We have traced this non-monotonous parametric dependence to the fleld dynamics 
taking place during the transition from the non-quadratic to the quadratic part of 
the potential. If the transition takes place soon after the beginning of oscillations, 
both a and a act as independent dynamical variables and, consequently, variation 
of the initial conditions can leave a non-trivial imprint on the expansion history. If 
the non-quadratic term always remains subdominant such features are absent. Also 
for the renormalizable case n = 0, the transition does not give rise to oscillatory 
behaviour but nevertheless leaves a clearly distinguishable mark in the results. This 
agrees with the analytical estimates we have derived for this special case. 

Some parts of the parameter space that give rise to an acceptable perturbation 
amplitude are hkely to be ruled out because they produce large non-Gaussianities. 
This is expected as we consider the subdominant limit r ^ 1 where the non-Gaussian 
effect are known to be large already for a quadratic model. However, because of the 
non-monotonic dependence C,{(^*) observed for non-renormalizable potentials, it is 
not evident what the flnal non-Gaussianity turns out be in this case. Indeed, when 
considering the expansion ( p.8|) for the curvature perturbation, the non-monotonic 
behaviour implies that the sign of (In a)" is varying, barring any accidental cancella- 
tions. In certain regions of the parameter space we should therefore flnd a signiflcant 
suppression for the non-linearity parameter /nl ~ (lna)"/(lna)'^ measuring the bis- 
pectrum amplitude. It is clear though that this qualitative argumentation cannot 
be pushed much further and a detailed survey of the parameter space is required 
to address the question of non-Gaussianity. This will be the topic of a forthcoming 
paper 
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